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Background The eikonal approximation is commonly used to calculate heavy-ion elastic scattering. However, the full evalua- 
tion has only been done (without the use of Monte Carlo techniques or additional approximations) for a — a scattering. 

Purpose Develop, improve and test the Monte Carlo eikonal method for elastic scattering over a wide range of nuclei, energies 
and angles. 

Method Monte Carlo evaluation is used to calculate heavy-ion elastic scattering including the center-of-mass correction and 
the Coulomb interaction. 

Results Angular distributions are presented for a number of nuclear pairs over a wide energy range using nucleon-nucleon 
scattering parameters taken from phase-shift analyses and densities from independent sources. A technique for the 
efficient expansion of the Glauber amplitude in partial waves is developed. 

Conclusion The center-of-mass and Coulomb corrections are essential. Angular distributions can be predicted only up to 
certain angles which vary with the nuclear pairs but all correspond to a momentum transfer near 1 fm"'^. 

PACS numbers: 25.70.Bc,25.55.Ci 



I. INTRODUCTION 



Monte Carlo (MC) methods are very valuable for per- 
forming multi-dimensional integrals, especially in quan- 
|;um applications [1]. Recently, calculations of nuclear 
structure using MC have been performed, both for light 
• • nuclei 0, y| and with the nuclear shell model [j. In 
. J^ the present article we study the use of this technique for 
K^ the evaluation of multi-dimensional integrals in the cal- 
^ culation of elastic scattering of heavy ions in the eikonal 
. . .approximation. 

Monte Carlo techniques have been used for classical 
simulations of nuclear reactions for many years (see e.g. 
Ref. [5| and references within). In that case (highly 
inelastic reactions), the assumption that phase informa- 
tion is unimportant appears to be valid, at least under 
certain conditions. This technique can reveal properties 
of nuclear reactions for inelastic reactions where many 
channels are open. Much less use has been made of the 
MC method for elastic (or nearly elastic) reactions. The 
reason for this might be that integrals performed with 
Monte Carlo are deemed not to be suitable to produce 
the highly oscillating amplitudes because of the growth 
in errors that typically occurs in this case. As we shall 
see shortly the evaluation of the Glauber elastic scatter- 
ing amplitude does not suffer from this problem. We now 



review some of the literature on this problem. 



A. Review 

1. The use of the Glauber approximation for 
nucleus-nucleus scattering 

Before discussing the use of Monte Carlo for eikonal 
calculations, we briefly mention the applications of the 
Glauber scattering method to heavy ion scattering more 
generally. 

Franco and Yin _6] and Yin et al. jj were the first 
to evaluate the full sum for the a — a scattering ampli- 
tude which appears to be the largest case for which the 
full scattering amplitude can be evaluated without fur- 
ther approximation and without the use of Monte Carlo. 
Their method relies on the use of a Gaussian approxi- 
mation for the nuclear density. Abu-Ibrahim et al. [8| 
evaluated the full Glauber multiple-scattering amplitude 
for proton-^He scattering. 

Al-Khalili, Thompson and Tostevin \9l studied the 
^^Li halo structure using eikonal methods, Al-Khalili and 
Tostevin 10] studied the relation of reaction cross section 
to radii in the Glauber model. They also studied proton 
helium halo nuclei citealkhalilil with Glauber scattering 
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and in Ref. [T2|, with Brooke, they treated -'^^B+^^C scat- 
tering with non-eikonal corrections, modifying the tra- 
jectory at lower energies. El-Gogary [l3| treated cluster 
nuclei to evaluate the Glauber formula. El-Gogary et 
al. [ij] used an approximate center-of-mass correction 
with a double Gaussian form for the density. Franco and 
Varma [l3 treated the center-of-mass effects to several 
orders and compared (primarily) with total cross sec- 
tions. El-Gogary et al. [la| did calculations with an 
approximate center-of-mass correction and noted that a 
simplified treatment of center-of-mass effects is problem- 
atic. Horiuchi et al. [l7[ did a systematic study of re- 
action cross sections usin g an improved expansion of the 
Glauber formula. Zhong 18] calculated nucleus-nucleus 
scattering in an a-cluster model with double Gaussian 
forms for the density. Charagi and Gupta [13| calcu- 
lated ^^O scattering from several nuclei using the Optical 
Limit Approximation (OLA). Abu- Ibrahim and Suzuki 
[20j used the OLA and a phenomenological profile func- 
tion to treat a large number of nuclei. Abu-Ibraham et 
al. [2l[ calculated halo nuclei scattering in the Glauber 
model. Abu-Ibrahim and Suzuki [22] calculated nucleus- 
nucleus scattering at intermediate energies using correc- 
tions to the OLA. Abu-Ibraham and Suzuki [2j] stud- 
ied the profile function in heavy ion scattering. Franco 
and Nutt [24] treated short-range correlations in heavy 
ion scattering. Lenzi et al. |2^, [2g| did a systematic 
treatment of heavy ion scattering in a OLA treatment of 
eikonal scattering. 



nuclei scattering with no approximation beyond those in- 
herent in the basic eikonal theory. No doubt this method 
should be used in any case that an exact calculation of 
the nuclear structure is available but a much more com- 
mon case is the scattering of heavier nuclei where such 
a solution is still for the future. In this case less ambi- 
tious approximations to the nuclear densities need to be 
used. We later will show calculations comparing the two 
methods and it would seem that the details (short-range 
correlations) make little difference. 



3. Center of Mass 

Of the above methods only that of Varga et al. [33] 
includes the center-of-mass correction exactly. We will 
see that this correction is very important, even for heavier 
nuclei. 

Chauhan and Khan [331 treated ^^C-^^C elastic scat- 
tering and found that center-of-mass effects play an im- 
portant role. The y u se an expansion of the profile func- 
tion. Liu et al. [3a] calculated a scattering from sev- 
eral nuclei at 1.37 GeV. They used an overall factor for 
the center-of-mass correction and varied the phase of the 
nucleon-nucleon interaction. 

Franco and Tekou [SJ] calculated the optical model ver- 
sion up to 5th order. They included corrections for the 
center of mass. Shukla [38] included the center-of-mass 
and Coulomb effects to investigate reaction cross sections. 



2. Monte Carlo 



4- Coulomb 



The Monte Carlo method has been used to evaluate 
the eikonal amplitude on several occasions. 

A. M. Zadorozhnyj, V. V. Uzhinsky and S. Y. Shmakov 
271 did perhaps the first application of eikonal Monte 
Carlo but no details are given. Merino, Novikov and 
Shabelski 28] very recently compared methods for ra- 
dius extraction. One of the methods is exact with Monte 
Carlo using the Metropolis algorithm [23] . Alkhazov and 
Lobodenka [30| calculated reaction cross sections for halo 
nuclei using the Monte Carlo method. Shmakov, Uzhin- 
skii and Zadorozhny [3l[ used Monte Carlo techniques for 
the generation of inelastic diagrams. Krpic and Shabel- 
ski [32j| used the Metropolis algorithm to calculate elas- 
tic and inelastic scattering from several nuclei using a 
diagrammatic expansion. They neglected the center-of- 
mass correction. Abu-Ibrahim and Suzuki [33] studied 
low-energy ^He-^^C scattering using the optical eikonal 
potential evaluated with Monte Carlo. 

Recently Varga et al. [SJj have combined the Monte 
Carlo Green's Function method of solution of few-body 
problems with Glauber theory to calculate a— a and halo- 



The basic Glauber method does not include the 
Coulomb interaction and a correction (often very impor- 
tant) must be made for it after the calculation. 

Kondratyuk and Kopeliovich i39] argued that it would 
be a good approximation to simply add the strong and 
Coulomb phases. Glauber and Matthiae [40] and Czyz, 
Lesniak and Wolek [41] used mainly this approximation. 
The corrections to this approximation were developed us- 
ing various techniques but most often based on an opti- 
cal potential. See Ref. [42] for a recent treatment and 
references to the previous work for the application to 
pion-nucleon scattering. Faldt and Pilkuhn '43\ used a 
semi-classical method with the Glauber model to correct 
for the Coulomb modification of the trajectories. This 
technique has the advantage that an optical model fit 
does not have to be made to the data to find the cor- 
rection factors. While they developed the method for 
pion-nucleus scattering it is currently being applied to 
heavy-ion elastic scattering. These two methods will be 



treated shortly is some detail. Charagi and Gupta [44 1 



treated low-energy heavy-ion scattering in the Coulomb 
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modified optical limit in two papers. Cha [45] introduced 
the deviation of the orbit in Glauber from the nuclear 
potential as well as the Coulomb potential. Charagi p6| 
commented on the Coulomb correction to the Glauber 
model introduced by Cha [43] saying that it failed to re- 
produce known reaction cross sections. Alvi et al. [43] 
calculated a— nucleus scattering treating the system by 
phenomenology (fitting to Ni) and showed the Coulomb 
effect. Ahmad et al. [1^] calculated ^^C-^^C elastic 
scattering using the first two terms of an expansion in 
a Coulomb modified Glauber calculation. 



B. Organization 

The paper is organized as follows. In section [II] the 
method is described showing how to satisfy the center-of- 
mass condition with the Monte Carlo sampling. In partic- 
ular, the technique for determining an auxiliary density 
which results in a desired center-of-mass single-particle 
density is developed. In this section the method for mak- 
ing the Coulomb correction is also presented. 

In section IIIII the calculations are compared with data 
for 20 nuclear pairs/energies. The parameters for the 
auxiliary functions are given. In section IIVI the results 
are summarized and conclusions are drawn. In the ap- 
pendices we present a method for the rapid partial-wave 
projection of the Glauber amplitude and explain how the 
variation calculation for the "^He density was done. 



II. MONTE CARLO EIKONAL 

NUCLEUS-NUCLEUS SCATTERING 

A. The Glauber Method 

The eikonal theory of elastic scattering as expressed by 
Glauber [49] and elaborated and studied by others J50H54| 
is believed to be an accurate representation of scattering 
at high energies and forward angles. It can be applied to 
calculate the multiple scattering of a simple projectile on 
a nucleus or collisions between nuclei, taking into account 
all of the scatterings possible in this theory. 

There are two major difficulties in applying this tech- 
nique: 1) as the atomic number of the scatterers increases 
the number of scatterings becomes very large and 2) the 
representation of the wave function (density) of the two 
scattering conglomerates raises the problem of the center 
of mass. The use of Monte Carlo techniques is well suited 
to solve both of these problems. 

For the scattering of a particles on ''He (the largest 
case treated to date for the full scattering series with- 
out using Monte Carlo) there are 16 possible scatterings. 
Taking into account all of the possible orders there are 



2^^ — 1 = 65535 terms in the multiple scattering expan- 
sion. They are not all independent however and, by di- 
viding them into classes, the calculation can be reduced 
[3] to 37 different types, each characterized by a 4 x 4 
matrix of zeros and ones to be calculated and included 
with different weights. There remains a 24 dimensional 
integral to be done for each term. By using a Gaussian 
representation of the wave function of "^He, Franco and 
Yin [6| and Yin et al. \J\ were able to provide expressions 
for these integrals and calculate the full sum. 

In the present paper the Glauber expression for the 
multiple scattering is used directly as a product of fac- 
tors in Eq. [1] without expanding into separate terms and 
the integral over the many-body nuclear density is done 
with Monte Carlo techniques, either with direct sampling 
or with the Metropolis algorithm. Normally one would 
be reluctant to use a Monte Carlo method to obtain a 
rapidly oscillating function as the amplitude for scatter- 
ing, but in this case Monte Carlo can be used to calculate 
the profile function (which is a smooth function since it 
consists of a sum of analytic functions) and then a stan- 
dard numerical technique can be applied to perform the 
last (one-dimensional) integral. The explicit development 
of the method is given in Section IIIBI 

The exact wave function for a nucleus at rest as a func- 
tion of the coordinates of the nucleons (we are ignoring 
spin for the present) would necessarily be translationally 
invariant in the absence of external interaction. In order 
to calculate the scattering between the two nuclei one 
would simply translate the wave function to place the 
center of mass at the origin. Of course one rarely has 
the exact wave function available but it is possible to in- 
troduce an approximate density which is, indeed, transla- 
tionally invariant where one can carry out this procedure. 
For "'He the four-body problem can be solved using Monte 
Carlo methods [3, Q • In the first instance a variational 
wave function can be used to find the minimum in en- 
ergy by varying the parameters in the trial function. The 
trial wave function should be translationally invariant. 
The square of the wave function (the density, which is all 
that is needed for the Glauber scattering calculation) is 
represented by a collection of Metropolis walkers. 

What has more commonly been done in practice is to 
start with an approximation to the density which is ob- 
tained empirically from a probe which is sensitive to the 
single-particle density relative to the center of mass, most 
often electron scattering. One then applies some trans- 
formation to assure that the center of mass of the nucleus 
is at the origin. Franco and Yin j6j] gave a prescription 
and showed that, using this formula, one could calculate 
with the fixed well assumption (for a harmonic oscillator 
potential) and then apply an exact correction. To illus- 
trate the importance of this correction we point out that 
the factor by which the amplitude is multiplied for a — a 
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scattering at — f = 3 (GeV/c)^ is 3.6 x lO'^. Clearly, a 
careful treatment of the center-of-mass correlation is im- 
portant. This problem is discussed at length in Section 
III El and the MC method is applied the evaluation of the 
Glauber amplitude for a — a scattering in Section IIII Al 
It is important to define what one means by a measure 
of the centcr-of-mass correction, that is, we need to de- 
fine the "zero effect" condition. It is common in both the 
double folding model and the OLA to take the density 
from electron scattering (corrected for the charge distri- 
bution of the proton) to form a product density. Often 
the center-of-mass requirement is ignored in these calcu- 
lations (see however Franco and Varma [l5| for a first 
order treatment). We will take this product density as 
our "no effect" model and compare the scattering from 
this case with scattering from a density with the proper 
center-of-mass properties. We will see that ignoring the 
center-of-mass effect, even for fairly heavy nuclei is ill 
advised. 



B. Basic technique 

The expressions for the Glauber amplitude are avail- 
able many places, see e.g. Franco and Yin [6j. 

The nuclear profile function, G{b) is the central func- 
tion in the theory and is given by 



5 = cr(l - ip)/(47ra) 

where the nucleon-nucleon amplitude has been approxi- 
mated by 



hi<l) 



ika{l-ip) „i, 



47r 



art. + lai, 



(5) 



where a is the nucleon-nucleon total cross section and p 
is the ratio of the real to imaginary part of the forward 
amplitude. For a discussion of the imaginary part of a, 
a/ see Ref. [50]. 

We have assumed a Gaussian approximation as a func- 
tion of momentum transfer for the nucleon-nucleon ampli- 
tude. The scattering parameters needed can be extracted 
from nucleon-nucleon (NN) scattering data [5g| and such 
a set is shown in Fig. [T] as a function of energy (see Refs 
|57| and [Sg] for tables of other determinations of these 
amplitudes. Of course, one might expect that these am- 
plitudes may well be modified in the nuclear medium (see 
e.g. Refs. [SO'-'GI'] but it may be useful to make calcula- 
tions with free values to see how large the corrections are 
likely to be. We take the free values used in the calcula- 
tions in this paper from the partial- wave parametrization 
of Arndt et al. [Sg. A plot of the free values used is 
shown in Fig. [T] For more details on the representation 
of the amplitudes, see the following section. 



„ A B 

G(b) = / \{ds,\{ds'^^\{s})i:\{s'})G{h,{s},{s'}), 

(1) 
In Eq. [l]the notation {s} represents the collection of all 

of the coordinates of the projectile or target and (/'({s}) 

denotes the projectile ground-state wave function and 

i/'ds}) the target ground-state wave function. 

A B 

G(b, {s}, {s'}) ^ 1 - n n [1 - ^'.(b + s. - s,)] . (2) 

The scattering amplitude is given by the two dimensional 
Fourier transform of G(b) and, since there is no direction 
defined for the nuclear wave function, it is just the Hankel 
transform of G(6), 

ik /■ f°° 

F(q) = — / cPhe'^^-^Gih) ^ik bdbJo{qb)G{b). 

(3) 
In Eq. [2] Tij (b -|- s^ — Sj ) denotes the two-dimensional 
Fourier transform of the elementary nucleon-nucleon am- 
plitude /(q), i.e.. 



rz,(b) = 2^Jd\e-^'i'^Mq) = ffe 



-\? jla 



(4) 



C. Nucleon-nucleon amplitude 

For NN free-space scattering there are essentially 4 in- 
coherent beams in the free (and unpolarized) case ac- 
cording to the 4 possible orientations of the spins of the 
two colliding nucleons. The imaginary parts of the for- 
ward elastic scattering amplitudes give the total cross 
section for each of these beams. The interference with 
the Coulomb in the forward direction can give the real 
parts of these amplitudes averaged with equal weighting. 

This combined amplitude is given by 



A 



1 



[M-t 



-++ 



M + A/+-+- + A/_+_ 



^M,.. 



M, 



= -[M, + M„] (6) 



with all amplitudes weighted equally. Here the ampli- 
tudes are labeled by "aligned" (subscript "a") and "op- 
posed" (subscript "o"). In the case of nucleus-nucleus 
scattering the amplitude would still consist of these two 
amplitudes if spin flip is neglected. Single spin flip am- 
plitudes are very small in the forward direction, so are 
often neglected. Since we are considering the scattering 
of two spin-zero nuclei, a single spin flip is not possible 



C Nucleon-nucleon amplitude 
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FIG. 1: Nucleon-nucleon parameters taken from the work of 
Arndt et al. [5o| . Also shown for the average NN cross section 
is the recent parametrization of ref. [6j| (dashed line). 



and this small amplitude must enter twice in the calcu- 
lation so that some other (correlated) nucleon can have 
its spin flipped in the opposite sense such that the total 
projection is again zero. These two constraints lead one 
to consider that the neglect of spin flip is a reasonable 
approximation. 

However, even in the forward direction there is an am- 
plitude M+ (. which does not vanish at zero degrees so 

a priori might be expected to contribute. However, we 
see that this amplitude must flip the spin of one nucleon 
in each of the nuclei. The resulting spin projection must 
be reduced to zero again by a second spin flip in each 
nucleus. This amplitude arises primarily from one-pion 
exchange. The average over spin (and isospin) removes 
this amplitude from consideration, at least in first order 
and we neglect it and consider only the amplitude coming 
from Eq. [B) 

The total cross sections which correspond to the two 
terms in Eq. |6] have been studied experimentally and 
their difference shows a rapid variation |62| sometimes 
attributed to a dibaryon (see e.g. (SJl)- If some mecha- 
nism could be found to change the equal weighting of the 



FIG. 2: Angular distribution of the nucleon-nucleon ampli- 
tudes taken from the phase-shift analysis of Arndt et al. [So ]. 
Also shown is the Gaussian approximation (dashed line) as 
used in Eq. [5] with parameters chosen to match the value and 
slope in the forward direction. The "aligned" amplitudes are 
shown in the upper panels and the "opposite" amplitudes in 
the lower panels. 



aligned and opposite amplitudes in Eq. [6] then a change in 
the energy dependence of the nucleon parameters could 
be expected. However, the spin projections come from 
two different nuclei so it is difficult to see how such a 
correlation might come about. The amplitude in Eq. |6] 
does not correspond to any directly measurable differen- 
tial cross section as a function of angle. Any attempt to 
extract the parameters directly from a measured differ- 
ential cross section will lead to parameters not suitable 
for calculations using the Glauber expressions; the am- 
plitude is a theoretical construct. Each of the terms in 
Eq. [6] can be well represented in the approximate form 
of Eq. [5] (see Fig. El. 

Listed in Table H] are values of the parameters used in 
the calculations presented in this paper. Also given are 
the maximum values of the center-of-mass angle and —t. 



D Nuclear profile function 



Case 


Ti,,6[Ref.] 


Ti„b/A (MeV) 


o" (mb) 


P 


aR 


ai 


Otra.ideg) 


-tMa. [{GeV/cf] 


«He-i2C 


0.230 GeVf64] 


38.3 


163.5 


0.680 


14.6 


-8.06 


20 


0.139 


^He-i^C 


0,250 GeV[65] 


41.6 


159.7 


0.737 


16.0 


-7.91 


10 


0.038 


a-'^O 


0.240 GeV 


[66] 


60 


103.2 


0.953 


20.8 


-6.65 


12 


0.051 


16o_16o 


1.120 GeV [67, 68] 


70 


86.4 


1.049 


23.0 


-5.87 


22 


1.224 


a-20«Pb 


0.288 GeV 


[69] 


72 


83.0 


1.069 


29.2 


-4.85 


33 


0.700 


a-^o«Pb 


0.340 GeV 


[69] 


85 


65.0 


1.169 


25.7 


-4.85 


30 


0.686 


12c_12c 


1.016 GeV 


[70] 


85 


65.0 


1.169 


25.7 


-4.85 


18 


0.560 


if^O-^Ca 


1.503 GeV 


m 


94 


60.5 


1.185 


26.3 


-4.56 


6 


0.254 


16q_12(. 


1.503 GeV 


[71] 


94 


60.5 


1.185 


26.3 


-4.56 


15 


0.565 


«-i«0 


0.400 GeV 


[72] 


100 


57.5 


1.196 


26.7 


-4.36 


37 


0.800 


a-20«Pb 


0.480 GeV 


[69] 


120 


47.4 


1.233 


28.1 


-3.70 


20 


0.443 


12c_12(^ 


1.449 GeV 


[73] 


120 


47.4 


1.233 


28.1 


-3.70 


12 


0.356 


12c_12c 


1.620 GeV 


[74] 


135 


43.1 


1.212 


28.6 


-3.84 


12 


0.398 


a-20«Pb 


0.699 GeV 


[69] 


175 


36.2 


1.095 


29.2 


-5.05 


12 


0.240 


12c_12c 


2.400 GeV 


[73] 


200 


34.3 


0.991 


28.7 


-6.42 


10 


0.410 


a«Ca 


1.37 GeV 


[75] 


343 


29.7 


0.435 


20.5 


-12.6 


12 


0.414 


a-i^C 


1.37 GeV 


[76] 


343 


29.7 


0.435 


20.6 


-12.6 


21 


0.799 


a- a 


2.554 GeV 


[73 


638 


38.9 


0.010 


10.5 


-8.10 


23 


0.765 


a-a 


4.20 GeV 


[78] 


1050 


42.5 


-0.303 


8.06 


-3.25 


13 


0.400 


a- C 


4.20 GeV [79, 8^ 


1050 


42.5 


-0.303 


8.06 


-3.25 


11 


0.716 



TABLE I: Parameters used in the calculations based on the nucleon-nucleon partial-wave amplitude analysis of Ref. (56 



D. Nuclear profile function 



Another advantage of the method is that the nuclear 
profile function (NPF), G{b), is available for study (see 
Fig. [3]). This function contains all of the information for 
the scattering since it is, essentially, the Hankel trans- 
form of the amplitude (we are dealing only with spin- 
zero on spin-zero scattering so there is only one ampli- 
tude). Clearly the NPF is strongly dependent on the 
single-particle density and the general shape reflects this. 

It is natural to ask how other characteristics of the 
nucleon distribution, such as short-, medium- and long- 
range correlations. They must manifest themselves in 
some manner but how? One can get some feeling for 
this effect by considering the "anatomy" of the calcula- 
tion. We see in Eq. [5] that the NPF is expressed as one 
minus a product. When this product goes to zero the 
NPF becomes unity which corresponds to total absorp- 
tion. When the impact parameter, b, is large only one 
factor (at most) in the product will differ from unity and 
only the single-particle density matters. As b becomes 
smaller, more factors differ from one and the product be- 
comes smaller in magnitude. Since the real part of the 
coefficient of the exponential in the function T, (7, is less 
than unity over most of the range of energies treated here 



each factor has modulus less than one. The imaginary 
part, a/, of a is smaller than the real part an and so 
plays a minor role in this qualitative discussion and we 
treat a as real for that reason. 

Just how small each factor is depends on the distance 
I Si — Sjl; if it is small then the factor is also small. This 
difference does not depend on nucleon-nucleon correla- 
tions since s^ and Sj are in different nuclei. However, 
since Real g is of the order of 0.25-0.8, no one factor can 
can drive the product to zero, it will take a combination 
of several factors. This will require that several nucle- 
ons be located close together and the probability of this 
occurring is sensitive to the correlations. Just how close 
together the nucleons have to be (in the 2-dimensional 
space) is governed by qr. Typical values of y^TiR are 
in the range of 1 fm. Repulsive correlations pushing the 
nucleons outside of this range would lead to greater trans- 
parency. Correlations shorter than this range can be ex- 
pected to have little effect. 

At the same time that the modulus of the product is 
decreasing it is developing a phase. Under the conditions 
just outlined the phase of each factor has the same sign so 
the phase grows monotonically as each factor is included. 
Thus the behavior of the NPF depends on the relative 
sizes of the phase and the modulus of the product. If 
the phase reaches n while the modulus of the product is 
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FIG, 3: The nuclear profile function for ^^C-^^C scattering 
with (solid) and without (dashed) center-of-mass correlations. 
The details of the calculation are discussed in section IIIIFI 



still sizable then the NPF will increase before becoming 
unity as seen in Fig, |3]with the CNC included. Further 
changes leading to the average modulus becoming larger 
(with only modification of the phase) could even lead to 
oscillations. 

In an effort to get a feeling for the effective number 
factors (thus the number of times that the correction is 
being applied) we calculated the profile function for val- 
ues of b in the crucial range of the surface as a function of 
the number of factors (see Fig, HJ , As we progress from 
large to small values of b more factors are important. For 
large b the profile function can be expressed as a sum of 
first-order terms and the dependence on the number of 
factors becomes linear. Of course in this limit there is no 
CMC since only first order in the density is involved. 
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FIG, 4: One minus the accumulated nuclear profile function 
as a function of the number of factors in the product in Eq, 
[2] for ^"Ca-^^Ca scattering. The labels on the curves are the 
values of fo in fm. 



center of mass which must be appropriately addressed. 
The Monte Carlo method can provide the solution. Con- 
sider two methods of treating the cm, motion. Each one 
provides a density in which the sum of all of the position 
vectors is zero but they are not (in general) equivalent. 
The first method relies on Gaussian densities and pro- 
vides a correction factor to a scattering calculation made 
with an auxiliary density centered about a fixed origin. 

Franco and Yin [61] used this first method, developed 
by Czyz and Maximon [8l| (see also Ref. [82|) which 
works only when using Gaussian densities because the 
center-of-mass can be expressed as a factor in this case. 
To implement this method, a model is made assuming a 
product of densities which are invariant under translation 
and using the algebraic identity: 



The scattering between two composite objects is 
treated in terms of the coordinate connecting their cen- 
ters of mass so that the particles making up the nuclei 
must be centered appropriately. This is the problem of 



g-Aa"fl"g-Q^(ri-R)"g-Q^(r2-R)"g-a"(r3-R) 



2 2 2 2 2 2 



Pa{ri)pair2)Pa{r3)--- (7) 
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where R = (^ ^i)/^ is the center-of-mass coordinate and 
Pa{r) is an auxihary density with reference to a fixed 
origin. The expectation vahie in Eq. [T] can be taken over 
the auxihary density (a calculation which is much easier) 
and the expectation value of the translationally invariant 
density obtained by dividing by the expectation value of 
the first factor on the left. 

For the density used by Franco and Yin with Pa{si) = 



IxiMr^ NS 






so that the auxiliary density has the same (Gaussian) 
form as the single-particle density but with a larger rms 
radius: r^ = Ar'^/{A — 1). 

We can also calculate the transform of the full density 
(Eq. [5]) in the general case as 



PsiM) 



dsids2 ■ ■ ■ dsA e 



N 



i(qi-Si+q2-S2---+qA-SA) 



dQY[pa{cir+Q/A) 



lx(M)l' 



N 
(2^ 






(8) 



where TV is a normalization factor. The single particle 
density relative to the center of mass, such as that ob- 
tained from electron scattering (see Refs. [SSl - ISa ] for cor- 
rections), for example, will be obtained, integrating over 
all but one of the coordinates, as 

Ps{si)= / rfs2ds3...|x({s})P 

= (^/rfQ e^^p4si)/-i(QM) (9) 

The Fourier transform of the single particle density (in 
the cm.) will be 



Psiq) = / dsi e^'i-^V.(si) 



N f 
j^ J dQ pa{q + Q/A)pt\Q/A) 



A^N 
J2 



N f 
^JdQpaici + Q)p^-\Q) (10) 



Normalizing, we have 

, , /rfQPa(q+Q)p^^(Q) .... 

J dQ p^(Q) 

For a Gaussian form in momentum space as Franco and 
Yin (y) assumed 

Paip) = e-P'/4"' (12) 

we can perform the integral on Q to find 

Ps{q) = e 4^«^ (13) 



_ /rfQntiPa(q» + Q) 

JdQp^iQ) 
Evaluating with a Gaussian form we have 



PsiM) 



-[E«?-i(Eq.)']/4o 



(14) 



(15) 



In the Monte Carlo method the algorithm leading to a 
density in the center-of-mass frame is to first select the 
A coordinates, u, according to A independent auxiliary 
densities, rjaiu), (assumed to be isotropic). The functions 
ria(u) have chosen forms and the principal aim of this 
section is to develop a method to pick the functions rja (u) 
such that they result in a specified single particle density 
relative to the center of mass. The coordinates to be used 
in the integration are obtained from the set of vectors u^ 
by 



lE"^ 



A 



(16) 



which means that the Monte Carlo densities that result 
from this transformation produce (by construction) func- 
tions which have the sum of their vector coordinates zero. 

To summarize, the procedure is to choose (isotropic) 
distributions according the density rjaiu) for all of the 
nucleons in a given nuclear configuration. The center-of- 
mass vector is then computed and subtracted from each of 
the Ui to give the coordinates to be used in the evaluation 
ofG(b,{s},{s'}). 

In order to choose an auxiliary function that gives a 
particular center-of-mass density it is very useful to have 
explicit expressions for the two functions that we need 
to connect. The Fourier transform of the single particle 
density relative to the center of mass will be given by 



/ 



Vsiq) — / duidu2 . . . duA X 
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^-1 \.A-l/'Z 



A-'^P'' VA 



(17) 



To compare with the resuhs of the previous section we 
can assume a Gaussian form again i.e. ?/a(p) — e~^ '^^ 
to get 



g^ (A-1)^ _ q^ A-1 A-1 q-' 



r]g(q) — e W^ a'^ e ^F'^^ 



'^ (18) 



which (taking a — /3) gives the same resuh as the method 
used by Franco and Yin. For the fuh Fourier transform 
we have 

r]s{{q}) = / duidu2 . . . duA x 



= / duidu2 ■ ■ . duA X 



■ ria[UA) 



e^Eli u..(q.-i Eti ^^)T^^{ui)r^^{u2) . 



If we take the Gaussian form again we have 



(19) 



(20) 



which, with a = (3 again, is the same resuh as found 
previously. So we see that for Gaussian functions the two 
methods give identical results. 

Eq. [T7]is general assuming all nucleons have the same 
density. For the case where the number of protons and 
neutrons is not equal and the neutrons and protons have 
different density shapes a different form holds. If Ca{q) 
is the Fourier transform of the proton auxiliary function 
and Cs{q) is the Fourier transform of the proton single 
particle in the center of mass (and r]a{q) and r]s{a) are 
the corresponding functions for the neutrons) then the 
functions will be related by 



Cs{q)^Ca[^<l)ctH^)V^{\) (21) 



Vsil) = Va {^^q) Vr\{)Cf {{) (22) 

If the auxiliary functions for the protons and neutrons 
have rms radii Rpa and Rna respectively the correspond- 
ing radii for the center-of-mass densities will be related 

by 

Rl, = {[{A- If + {Z~ l)]Rl, + NRl,} /A' (23) 



i?L ^{U- If + {N- l)]i?L + ZRl,} /A' (24) 

For the case of a given form chosen for the single- 
particle density, we need to find an auxiliary density 
which satisfies Eq. [TTlor 



Vaiq) Va 



A~\ 



A-\ 



= Vs 



_Aq_ 

A-1 



(25) 



One can find a solution to this equation for small q, by 
first expanding both ?/s(g) (assumed to be known) and 
r]a{q) for small q. 



?7a(g) = 1 -/V---; Vs{q) = 'i^ - i^q'^ ■ ■ 

Then we have to first order in q^ 

(l-f,q^)[l-f,qyiA-in^-'^ 

= (i-w2)[i-mV(A-i)] 






(26) 



so that 



^J■ 



vA 
A~\ 



(27) 



(28) 



and the same relation holds between the rms radii in the 
general (Monte Carlo) case as for the Gaussian case. 

We could (in principle) solve Eq. [55] numerically on a 
mesh by constructing the first few elements on the mesh 
(3 of them) with the expansion. Then we can continue to 
calculate the remaining values on the mesh by evaluating 
the \th point on the mesh with 



■^aisii 



Vs 



_Aqi_ 
A-1 



l< 



A-1 



A-1 



(29) 



where the value of r]a in the denominator is obtained 
by interpolation from previously calculated values on the 
mesh. These values would always be available since the 
argument is much smaller. This process works well for a 
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Gaussian form but is not stable when there is a zero in 
the form factor. 

In practice, we choose parametrized forms for rjaiQ) 
and vary the parameters using Eq. [T7] to fit an (assumed 
known) function, ris{q). For convenience we choose func- 
tions to represent ria(r) which can be readily sampled 
directly (see e.g. Ref. (87|). 

The density with the center-of-mass effect is distin- 
guished from the product density with the same single 
particle radial distribution by the fact that the distances 
between members of any pair of particles is greater than 
it would be for a simple product density. We can see this 
as follows. 

For a single-particle density with independent particles 
the average of square of the distance between particles is 
given by 

< (ri - r2f >= 2 < r^ > -2 < ri • ra >== 2 < r^ > 

since the average over the independent particles gives a 
cross term of zero. If we include a center-of-mass condi- 
tion 



Er.-O 



(31) 



i=l 



then, replacing r2 in the cross product 



< ri •r2 > = 



< r' 



>-E<'"i 



r,: > 



<r^ > -{A -2) < ri -ra >= 



Thus, with center-of-mass correlations 



< (ri -r2)^ >=: 2 < r^ > 



< r^ > 
A-1 



A 



A-1 
For individual nuclear configurations we have 



E 



»#j^ 



r.f 



2A E 



k ' k 



A{A-1) 



1 



(32) 



(33) 



(34) 



Thus we see that there is a correlation among pair of 
nucleons arising from the cm. corrections. This condi- 
tion is realized explicitly in the present method since the 
auxiliary density has a squared radius which is a factor 
of A/ {A — 1) larger than the single particle density in 
the center of mass and, since the positions are indepen- 
dently chosen, Eq. [30]holds with the larger square radius. 
Since the displacement of the density to satisfy the center- 
of-mass condition does not change the distance between 
nucleons, this larger inter-particle distance is preserved 
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FIG. 5: The ratio of inner-|-outer corrections to outer cor- 
rections only for the Faldt-Pilkuhn [43] (point Coulomb) and 
optical model with a realistic charge density. The full angular 
distributions are given in section UlI HI 



while the distance of individual particles from the center 
of the nucleus is reduced. 

With the procedure used here there is an initial den- 
sity constructed in a Monte Carlo sense and then each 
realization is shifted by an amount to put the center of 
mass at the origin. In a spherical density with each par- 
ticle thrown independently of the others, the rms dis- 
tance between any pair of particles is \/2 times the rms 
radius of the nucleus. Since the shift of the entire nu- 
cleus does not change the distance between pairs and the 
initial (auxiliary) density has a larger extent than the 
final density relative to the cm. then the relative dis- 
tance between pairs in the final density will be larger 
than the one which would be associated with a density 
constructed from independently thrown nucleons with the 
shape of the center-of-mass density. As the mass num- 
ber, A, goes up this effect will become smaller with the 
relationship for the inter-particle radius squared being 
RcMC — ^/(^~ 1) ^ ^No CMC- Hence one might assume 
that the center-of-mass effect goes to zero as A increases. 

While it is true that the basic effect is becoming 
smaller, the calculation of a given observable may not. 
In the case that we are treating, of multiple scattering, 
the number of scatterings increases with A. Each scatter- 
ing among the nucleons depends on the relative distance 
between nucleons and, while this distance is approach- 
ing that which would come from an independent particle 
density, the smaller effect is applied more times (there 
are more scatterings) so it is a numerical question as to 
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whether the effect decreases or perhaps even increases 
with A. If the fundamental scattering interaction is weak 
so that only the low order scatterings are important then 
one can expect that the effect decreases with A. How- 
ever, for the strong NN interaction it is possible that the 
importance of the CMC to multiple scattering does not 
decrease with A at all. Consider, for instance, the case of 
the scattering of a nucleus with A nucleons on its twin. In 
this case the number of scatterings goes as A^ although 
there will be other factors which will limit the effective 
number of scatterings. 



F. Sampling considerations 

The multidimensional integral in Eq. [l]is to be done by 
Monte Carlo. The function to be averaged over, [G(b)], 
is complex but, as we shall see, the densities can be easily 
sampled. Ordinarily Monte Carlo is not very accurate if 
one takes a Fourier transform over a binned distribution 
(which we are not doing here). The integral in Eq. [3] is 
done by a standard quadrature method. 

It is convenient to sample from pools of configurations 
of nuclei. For a variational or Monte Carlo Green's Func- 
tion method for generating the nuclei this is the natural 
way to carry out the sampling process, but even if direct 
sampling is being done it is a useful method. One creates 
Nf configurations for the target and Np for the projectile. 
Then M samples are taken in pairs, one from each pool, 
of all of the nucleons for a nucleus from the same configu- 
ration. If M is large enough the pairs will be identical for 
some cases. However, the number of identical pairs will 
be M'^/NtNp for a fraction of repeated pairs of M/NtNp. 
Thus for a typical case of one million entries in each pool 
and 10 million Monte Carlo samples the pairs repeat only 
10'^/10i2 ^ 1/10^ fraction of the time. 

This technique is also useful for calculating the effect 
of no center-of-mass correlations (CMC). One can simply 
take the coordinates for each nucleon from a different con- 
figuration. This guarantees that the same single-particle 
center-of-mass density is used but the coordinates in each 
nuclear configuration are uncorrelated. 

The auxiliary functions are chosen with a form which 
can be easily sampled. A common form used in many of 
the cases which follow is 



by 



. , un+l - 

£ — ^ rj 1 



(35) 



Often only two terms in the sum are needed to get an ade- 
quate representation. To sample this function one term is 
selected with probability, (7„, and then the corresponding 
normalized probability function is sampled. The func- 
tions in this case can be sampled easily by choosing r 



innii^i^x 



(36) 



where the set, AT^, are independent random numbers uni- 
formly distributed between and 1. 



G. Coulomb correction 

If the full amplitude with Coulomb is 

-J oo 

^(^) = VI E(2^ + 1)(^^ - l)^^(cos0) 



^-0 



. oo 

= !r{Q) + ^^(2£+ 1)(5^ - e2'"0A(cos^ 



£=0 



-. oo 

/fl,(e) + T-T ^(2£+ l)e2-^(5,e-2-. _ i)P,(cos( 



£=0 



(37) 

where jn(&) is the Rutherford amplitude and the ag are 
the Coulomb phase shifts, and the amplitude without 
Coulomb is written 



_. oc 



(38) 



£=0 



then we can write the factor between the two S-matrix 
elements as 



Co — S, 



ee 



-2iu 



'isl 



(39) 



(see Ref. [42 and references therein). The inner Coulomb 
correction can be calculated with the use of an optical 
model. The procedure is to calculate Ci from Eq. [39] 
using for Sg the result of an optical model fit with full 
Coulomb and for 5^, the value from the same optical 
model without the Coulomb interaction. The quantity 
Ci is then used as a correction factor with S^ coming 
from a partial wave expansion of the Glauber amplitude 
and Sg being the Coulomb-corrected result. 

Another way of getting the inner correction was given 
by Faldt and Pilkuhn ^4^] as a simple shift in the value 
used in the integral over the profile function. 



where 



F(q) = ik / bdbJo{qb)G{b') 
Jo 



b' = ^JW^^rff^ + 77/ A:, 



(40) 



(41) 
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77 is the Coulomb parameter rj — ZZ'ac/v and k is the 
center-of-mass momentum. The Coulomb corrected S- 
matrix element will be given by 



Se 



2iS, 



SFP^+2^., 



(42) 



where 5*^^ is the matrix element coming from the ex- 
pansion of Eq. |3D]in partial waves. For a similar method 
of treating the Coulomb correction see Vitturi and Zardi 



In order to carry out either of these corrections it is 
necessary to have the amplitude expressed as a partial 
wave sum. Since the eikonal calculation gives the am- 
plitude as a function of momentum transfer or angle the 
projection of the amplitude onto partial waves is needed. 
The projection was not made directly but by using Eq. [3] 
we can write 
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FIG. 6: Ratio of calculations for a — a scattering for different 
numbers of Monte Carlo evaluations showing the degree of 
calculational precision obtained. The a-particle density used 
is the same one used by Franco and Yin [g] 



Se-l = -k^ dx bdbPe{x)Jo[q{x)b]G{b). (43) 
J-l Jo 

The integral over x of Pe{x)Jo [q{x)b] can be carried out 
first very efficiently (see appendix A). The result was then 
integrated over b with b G{b). As a check the amplitude 
was then reconstructed from the S-matrix elements and 
compared with the original calculation. With the method 
given in Appendix A agreement was found within 0.1%. 

For the case of ^^O-^^O scattering at 1.12 GeV we com- 
pared the two methods. Figure [5] compares the two meth- 
ods of making the inner correction for ^^0-1-^^0 (see Sec. 
IIII HI for details). It displays the ratio of the fully cor- 
rected amplitude to that with only the outer correction. 
It is remarkable how well the two methods agree. The 
disagreement beyond 16 degrees is due mainly to the fact 
that the optical model is calculated with a finite charge 
distribution and the method of Faldt and Pilkuhn as- 
sumes a point charge distribution. For other treatments 
of the Coulomb correction, see Refs ^] and jMj . 



III. APPLICATION OF THE MONTE CARLO 
METHOD: RESULTS 



In this section we carry out calculations using the full 
Glauber formalism with values of the parameters taken 
from fits to the nucleon-nucleon amplitudes [5a| ■ In some 
cases we will vary those parameters to investigate the 
sensitivity to their values. Table U summarizes the cases 
treated, the data sources and the free-space parameters 
used. 



A. a — a scattering 

There are some interesting points for the a — a scatter- 
ing calculation. Notice that there is no separate center-of- 
mass correction factor as there was in the formulation of 
Franco and Yin [6j . The Monte Carlo method calculates 
directly with the many-body density with the method 
outlined above. As a test of the method we repeat the 
calculation of Franco and Yin [6| using their technique 
and the Monte Carlo method using Gaussian densities. 
From the previous discussion on the center-of-mass cor- 
rection we should expect to find the same result aside 
from Monte Carlo statistical errors. Figure [5] shows the 
ratios of the calculations done in the two ways. It is seen 
that the results are the same to about 2% even though 
the absolute magnitude of the cross section changes by 
almost 10 orders of magnitude over the angular range 
considered. The difference between the two calculations 
is invisible on a log scale. For 10® Monte Carlo evalu- 
ations the calculation took 10 hours on a standard 2.5 
GHz PC. 

The nucleon density of ^He has been somewhat of 
a puzzle. The fits to electron scattering [91| with the 
method of a sum of Gaussians (SOG) show a large de- 
pression in the center of the nucleus that the Monte Carlo 
Green's Function calculations do not find 3]. However, 
the x^ per data point is considerably less than unity show- 
ing that perhaps an over-fit was present in the SOG fit. 
Using other forms to represent the density, acceptable fits 
with a significantly smaller depression can be found. We 
have made a fit to the electron-scattering data [9l| with 
a density in which we constrained the ratio of the density 
at the origin to that at the peak to be 0.4 in order to have 
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Points: Pudliner et al. 
Solid: Variational density 
Dash: Fit to MSW 0.4 
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FIG. 7: Comparison of three densities used for the a particle 
in the calculations. The reference MSW is to McCarthy, Sick 
and Whitney |91|. The points were read from Fig. 15 in 
Pudliner et al. 



a density similar to that found in Ref. [9l| but with a 
less severe depression. That density is shown in Fig. [7] as 
the dashed line. We fit this density to find an auxiliary 
density given by 

r^PMSwir) = (1 - a)N, {e'^^^ - 6-"='')'° 



+ a[0.95N2e-'''-''"^'^''« + OmNse''-'-^'''^'^'''] (44) 

where a = 0.43 and 

Ce = O.Slfm^ ; be = 0.61fni^ ; 

ra = 1.28 fm; n = 1.8 fm; 



oo = 0.1 fm; ai = 0.3 fm 



(45) 



Ni , N2 and N3 are chosen to normalize each of the indi- 
vidualprobability densities. The calculations of Pudliner 
et al. [3| give perhaps the best estimate of the ^Hc density 
and show only a slight depression in the central region. 
Also shown in FigJ7]is a fit to the Monte Carlo density 
of Pudliner et al. [3] (large dots) from reading the points 
from the graph. The auxiliary density for this fit is given 

by 

r^Pud = X^d^^r^^e-'^-'/Ul + (1 - Xi)fr^e-f'-/8\ (46) 

where A4 = 0.835 and the resulting density is shown in 
Fig. [7] (dotted line). In order to have a density with short 
range correlations included we have carried out our own 
variational calculation (See Appendix B). 



FIG. 8; Calculations of a — a scattering showing the effect 
the center-of-mass correlations. The dash-dot curve uses the 
same Metropolis density as the solid curve but each nucleon 
is drawn from a different (random) realization of the nucleus. 
Thus it has the same single particle density but the effect of 
the CMC is missing. The data are from Berger et al. [77| 



To test for the relative importance of the short-range 
correlations and the center-of-mass correlations, we per- 
formed the calculation as outlined in the previous section 
starting from the single-particle density resulting from a 
density derived from a binning of the variational result. 
The cm. density from that calculation is also shown in 
Fig. [7] (solid line). 
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FIG. 9: a — a scattering at TLab ~ 4.2 GeV showing the 
dependence on the density used. The notation MSW refers 
to the fit to the electron scattering fit by McCarty, Sick and 
Whitney [9l|] as described in the text. The partial wave anal- 
ysis (PWA) is that of Ref. [5a|. The data are from Satta et 
al. ;78] 

The auxiliary density, r]a{r), as fit to the variational 
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density can be expressed as 



r'^m&r{r) == (1 - a)gi{r) + ag2{r) 



(47) 



where 



51 W 



Ngr^e-^o' 



-bcr\2 



32 = N.ie-'^'-' ~ e-''-) 



a = 0.505, ag = 0.56 fni~% 
= 0.779 fm~\ he = 1.345 fm"^ 



(48) 



where Ng and TVe were chosen to normahze gi and 52 in- 
dividually. Note that the volume element, r^, is included 
in each case. Since the wave function has an asymptotic 
limit of the form e~'""/r then, with the volume element 
included the tail should behave as an exponential. 
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FIG. 10: a — a scattering at 4.2 GeV showing the variation 
with the parameter p. The data are from Satta et al. [Ta ]. 
PWA means the partial wave analysis from Ref. [5g] 

Figure |8] shows a test of the importance of short-range 
correlations and center-of-mass correlations. For this il- 
lustration only, the parameters have been chosen to give a 
reasonable representation of the data, unlike most other 
figures presented which use the free parameters deter- 
mined from amplitude analyses [56|. The solid curve 
shows the calculation with the full variational wave func- 
tion and the dashed curve shows the results obtained with 
the auxiliary density used to correct for the center of 
mass. The dash-dot curve shows the result of choosing 
each nucleon from a different configuration of the nucleus. 
In this case the single particle density will be identical to 
the other two cases but there is no CMC among the nucle- 
ons. The principal difference between the two is that the 
solid curve does not have the short-range correlations that 
are in the variational wave function. Hence, at least for 
the variational calculation performed here (see Appendix 
B), they are not important in agreement with UUo and 
Feshbach (H. 
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FIG. 11: The top panel shows the fit to the C charge density 
using the auxifiary density in Eq. 1491 The electron-scattering 
data are from Ref [93]. lower panel shows the point density 
along with the auxiliary density. Also shown in the lower 
panel are the points obtained from binning the radius values 
in the Monte Carlo calculation, the open circles corresponding 
to the auxiliary density and the sofid circle to the center-of- 
mass density. 



Figure [HI shows the dependence on the density used. 
It is seen that the sensitivity is very small and reason- 
able variations would not seem to be able to significantly 
improve the agreement with the data. 

Figure [TOl shows a comparison of the calculation at 4.2 
GeV beam energy with the data of Satta et al. [78| . Also 
shown are the results of a variation with the p parameter. 
This parameter is the least well determined experimen- 
tally of the nucleon-nucleon parameters. It is seen that 
the prediction is rather poor when compared with the 
data and that moderate variations of p are unlikely to 
improve the agreement. 



B. Q-^^C scattering 

The carbon auxiliary density was obtained from the 
modified harmonic oscillator fit of charge density ob- 
tained from electron scattering data [92| . The form found 
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taken from Pudliner et al. [3| as described in the text. 

for the auxiliary density is 

r^r]i2ir) = 0.90916 hi{r) + 0.09084 /i2(r) (49) 

where 
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hi{r) = Mir'-^e^'"''"'; /i2(r) ^ Ahie-"^^ 



-bi2r\2 



(50) 

and di2 = 0.7311 fm-\ aia = 1.7488 fm-i; 612 = 0.6429 
fm~^. Ml and M2 are chosen to normalize hi and /i2 
separately. 

Figure [11] shows the auxiliary density and the center- 
of-mass body and charge densities which result. For the 
light nuclei the auxiliary density is normally quite differ- 
ent from the center-of-mass density. 
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[73 ] showing the effect of the Coulomb correction. 



Figure [T2] shows the variation of the differential cross 
section at 4.2 GeV with the parameter p. As for a-a 
scattering at the same energy (Fig. [9]) the moderate vari- 
ations of p do not really improve the agreement with the 
data. Figure [T^j shows the prediction of the calculation 
for a-^^C scattering at 1.37 GeV kinetic energy with the 
parameter a/. The agreement is no better than that at 4.2 
GeV and variations of the other parameters give similar 
results. Again, it seems unlikely that a modest variation 
in parameters will bring the calculation in line with the 
data. 



C. a- O scattering 

The auxiliary function for ^^O was taken to have the 
form: 

r'^Tjieir) = 0.8646 vi{r) + 0.1354 V2ir) (51) 

with 

vi{r) = A'ly^e-'^''"'; V2ir) = Nhie-"''"' - e'^'^^'f 

(52) 
and 

ai6 = 2.738 fm-^ bie = 0.2976 fm"^; dig = 4.174 fm-^ 

(53) 
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FIG. 15: 0-^*^0 scattering at 240 MeV. The data are from 
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FIG. 16: Q-^^C scatter ing at 1.37 GeV compared with the 
data of Alkhazov et al. [7511 



with 



and 



Ml = d}J/10! 



Ah = 2ai66i6(ai6 + bi6)/{ai6 - he)^ 

chosen to normalize vi and V2 ■ 

Figure [T3] shows the results of the calculation a-^^0 
compared with the data of Wakasa et al. ff^ at 400 



MeV. Here the agreement at forward angles is satisfac- 
tory. Figure [15] compares the calculation at 240 MeV 
with the data of Lui et al. [66] here the agreement is 
fairly good (aside from a slight normalization problem) 
up to about 7 degrees. After that there is a consider- 
able difference. The Coulomb correction is substantial, 
especially at larger angles. 



D. a- Ca scattering 

The auxiliary density for Ca was obtained from a point 
density extracted from pion scattering [9J| which is very 
similar to that obtained from electron scattering |95| . 
The form is 



r^Tjioir) = 0.740072 vi{r) + 0.259928 waW 



(54) 



with 



vi{r) = Mir^^e-'^^'^; V2{r) = MzCe""'""' - e-'"""^)^ 

(55) 
and 



a4o = 0.641 fm" 



t^40 



0.443 fm" 



^40 = 3.597 fm~' 
(56) 
with Ml and M2 chosen to normalize vi and V2- 

The results of a calculation for a-'^^Ca scattering at 
1.37 GeV are shown in Fig. [TSIand compared to the data 
of Alkhazov et al. [73. The Coulomb correction plays 
a very large role. There would be no agreement at all 
without it. While the agreement is far from perfect, it is 
satisfactory for no adjustable parameters at least for the 
angles less than 8 degrees (with the exception of the two 
forward data point which have large errors). For other 
treatments of a— ''"Ca scattering see Refs. [9g| and [93] • 



E. a-^°*Pb Scattering 

The auxiliary density for ^"*Pb was fit to the charge 
density of ^°*Pb so corresponds to the proton density. 
It is believed that the neutron density is different from 
the proton density but a study of the effect of a different 
neutron density, interesting though it might be, is beyond 
the scope of the present work. Since the lead nucleus is 
much larger than the other nuclei considered so far, with 
a significant flat portion inside the surface, the forms of 
the auxiliary density used for the lighter nuclei are not 
appropriate for this case. Because of the large number 
of nucleons, the difference between the auxiliary density 
and the center-of-mass density is not very great. For 
these reasons a Woods-Saxon (WS) density was used, fit 
to the charge density [93]. The volume element must 
be included as well in the sampling as has been done in 
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FIG. 17: Q- Pb scattering compared with the data of Bonin et al. [69 



previous cases. While the WS density can be directly 
sampled, it is the form 



(57) 



l + e(^) 



which is needed. This sarnpling was done by the method 
of selection of variables [871 ]. A comparison was made 
between variables from two samples. First, ri is chosen 
according to a linear distribution from to Rq then a 
second value, r2 was obtained from the WS distribution 
with 



7'2 = aln 



1 



b[{i + i/by 



where -F is a random number uniformly distributed be- 
tween and 1 and b — e~'^^°- . If r2 is greater than r\ it 
returned as the desired value of r. Otherwise the process 
is repeated. The aparticle density used was the one fit 
to Pudliner et al. [3 as discussed earlier. 



Figure [17] shows the results of calculations at energies 
of 288, 340, 480 and 699 MeV compared with the data of 
Bonin et al. [63] • The Coulomb correction plays a very 
large role in determining the cross section as might be 
expected. The agreement is rather good in the forward 
direction but worsens rapidly beyond a certain angle, dif- 
ferent for each energy. 
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FIG. 18: ^^C-^^C scattering at 1.016 GeV with and without 
the center-of-mass correction. The data are from Buenerd et 
al. Izi 



F. i^c.iac Scattering 

A comparison with the data at 1.016 GeV [70] and 
1.45 and 2.4 Ge V i73il has been been ma de by several 
groups P.[ll.[20l.[2if2l.[25l.[2l.l4g.[99l4l0TI| using various 
met hods . See also the phase-shift analysis by Mermaz et 
al. |l02l |. We can use the -^^C densities found in section 
1111 Bl to calculate the scattering of ^^C from ^^C. 

We have carried out a test of the size of the center-of- 
mass effect in much the same way as was done for the 
variational wave function in a — a scattering earlier. The 
^^C-^^C scattering calculations just presented were made 
by repeating the realization (in a Monte Carlo loop) of 
two carbon nuclei and then carrying out the evaluation 
of the necessary equations. The calculation can also be 
done by first constructing a pool of nuclei, each properly 
centered about the center of mass (one million were used 
in the current calculation) and then drawing complete 
nuclei randomly from this pool for each carbon nucleus. 
These two methods give the same result. 

One can now modify the calculation, in the same man- 
ner as before, so as to choose the 12 different vector co- 
ordinates for each carbon nucleus from 12 different re- 
alizations in the pool. In this way one is guaranteed to 
have the same single particle density but with uncorre- 
lated particles. The result of a calculation for this energy 
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with a spherical density for ^^C is shown in Fig. [TSl where 
it is seen that both the Coulomb and center-of-mass cor- 
rections play a large role. The agreement is good up to 
an angle of about 5 degrees {t « —0.045 (GeV/c)^) but 
poorer after that. 

The result without center-of-mass corrections is shown 
as the dashed curve. While this curve is shown as the 
ratio to Rutherford as a function of angle and Fig. [5] for 
helium gives the absolute cross section as a function of 
t, one can see that the effect is very similar in the two 
cases. Thus, at least for this very limited sample of two 
cases, the center-of-mass effect does not decrease. 

Figure 13] shows the profile function for ^^C for the cases 
with and without CMC. One can understand the behav- 
ior to some extent. For small values of b the many factors 
in the product in Eq. [U (144 for the case of ^^C-^^C scat- 
tering) , a large fraction of which have magnitude less that 
unity, will cause it to be very small. This very small cor- 
rection to unity leaves the profile function at essentially 
one in this region of b, hence is insensitive to the CMC. 

On the other hand, for large b one can expand the 
product in terms of single, double, triple etc. scatterings. 
Since a double scattering will have two factors of the 
Gaussian function, the triple scattering will have three 
factors etc., the single scattering will dominate for large 
b. Since we are holding the single particle density in the 
cm. fixed, the CMC will have no effect in this region. 
Because the large b values dominate the forward scatter- 
ing we must expect the small angle cross section to be 
insensitive to the CMC as is observed. 

Thus, it is only in a relatively small region of values of 
b that the effect will be influential. From Fig. [3] this is 
for 3.5 < b < 5.5 fm. Since, as one decreases b from the 
external region the sensitivity to the CMC will become 
greater as the number of scatterings goes up until this 
number gets to be so large that the product becomes 
very small in magnitude. 

It is generally believed that carb on ha s a strong oblate 
deformation. Lesniak and Lesniak [l03| included this ef- 
fect in proton-carbon scattering using Glauber theory. 
While they developed the proper, fully quantum, theory, 
in the end they used a semi-classical approximation sim- 
ply avera ging the amplitude over rotated densities (see 
also Ref. Il04l |). Some recent studies (of fission and reac- 
tions [1055) have treated the problem in the same way. As 
pointed out in Ref. [103], the electron scattering (at least 
in the single interaction approximation) should not de- 
pend on the deformation so the average density remains 
the same as before. We assume a form for the density 
symmetric about the z-axis and with a distribution in 
the polar angle given by 



p{r,6) ex po{r)ain"^ 9. 



(58) 
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FIG. 19: Comparison of three deformations of the carbon nucleus at four energies. The dashed curves correspond to a spherical 
nucleus (n=0), the solid curve to a ratio of 1.14 (N=2) and the dash-dot curve to a ratio of 1.73 (n=4). The data are from 
Buenerd et al. JQ,], Hostachy et al. [73] and Ichihara et al. [74| . 



Thi s sim ple form is inspired by considerations from 
Ref. Il06ll fpage 62). The paper of Svenne and Mack- 
intosh [lOTJ I presented arguments why ^^C was known 
to be deform ed in response to the paper by Friar and 
Negele llOS^ who pointed out that, with the usual form 
of the deformed density, the existence of a deformation 
for ^^C was contrary to the electron scattering measure- 
ments since the fall off of the density in the surface region 
was strongly affected. The form of Eq. [SUdoes not suffer 
from this problem as can be seen by taking the example 
of n = 2. In this case we have 



Po{r) 



^ _ ^ V r'^ 

3 15 -^ ^ 



K)Yn^a,^a)\ (59) 



where the Legendre polynomial has been expanded in 
terms of angles relative to a fixed axis. Here {9p,(l>p) 
are angles of a given nucleon in the nucleus with respect 
to a fixed axis and the angles {0a,4'a) are those of the 
body symmetry axis relative to the same fixed axis. With 
a one-body operator the second term will give a ^2(9^) 
transform of the density but when the average over the 
direction of the body axis of the carbon nucleus is taken 
it will vanish. Hence a one-body operator probes only 
the "spherical part" of the density. 

For n = 2 the probability in 9 is given by 



p{r,9) ex po{r)sin^9. ex -poir)[l ~ P2(cos( 



PiO) 



■ sm 



(60) 
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< x^ > / < z^ >== 2 
For n = 4 the probability in 9 is given by 
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(62) 



9 1 o 9 9 O 9 
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< a;^ > / < z^ >= 3 



(63) 



These densities look more like donuts than ellipsoids 
of revolution (assuming that this density is applied to all 
the nuclcons and not just the p shell). We note that a 
cluster models would also lead to densities with zero at 
the center. We can estimate the /32 parameter describing 
deformation by using prescriptions given by Hagino [1091] 
based on the rms radii. 

We take 



1 



\P^^ 



< Z'^ >2 



SO that 




(64) 



(65) 



For n = 2 we have P2 ~ —0.34 and for n = 4 
j32 — —0.52. Svenne and Mackintosh [lOTJ gave a sur- 
vey of values obtained for /?2. From their list we see that 
the Nilson model gives P2 = —0.64 and the a c luster 
model /3 — —0.41. A little later Vermeer et al. [llO| mea- 
sured the quadrupole moment of the 2+ state from which 
one can infer a quadrupole moment for the ground state 
of —22 ± 10 e fm^ implying a value of /32 = —0.57 with 
a 50 % error. In a recent calculation of ^^C-^^C fusion 
(l05| the value /32 = —0.4 was used (along with a hexade- 
capole term) so these forms lead to reasonable values of 
the deformation. 

The results for n=2 and n=4 are shown in Fig. [TOland 
compared with data. Again it is seen that agreement 
with data is good in the forward direction but deterio- 
rates rapidly at larger angles. The deformation plays a 
significant role at larger angles. 
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FIG. 20: The ^''O-^^C scattering cross section showing 
effect of the Coulomb and center-of-mass corrections, 
data are from Roussel-Chomaz et al. uM 



the 
The 



G. if^O-^^C Scattering 

The scattering of Oxygen on Carbon is shown in Fig. 
[201 The effect of the correction for the Coulomb inter- 
action is quite large. For the same scattering the effect 
of the correction for the center of mass is also shown. It 
is seen that the agreement a very forward angles is very 
good with the Coulomb correction playing a substantial 
role. While it cannot be said that the agreement with 
data is good at larger angles, the calculation follows the 
trend of the data quite well which is not true if either the 
Coulomb or center-of-mass correction is not included. 



H. 



^0-"*0 Scattering 



Using the density found earlier we now calculate Oxy- 
gen scattering from Oxygen at 1.120 GeV. The compar- 
ison with the data of Nuoffer et al. [63| is shown in Fig. 
[5TJ For angles beyond 5 degrees there is a considerable 
discrepancy (more after the Coulomb correction) but in- 
side of that angle the agreement is quite good. 

In order to test the inner Coulomb correction we made 
an optical- model fit to the data [67|. A simple form was 
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FIG. 21: I'^O-^^O scattering at 1.12 GeV. The data are 
Nuoffer et al. [13, and Khoa et al. B 
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FIG. 22: The ^®0-*°Ca scattering cross section showing the 
effect of the Coulomb and cm. corrections. The data are 
from Roussel-Chomaz et al. uM 



used 



Kpt(r) 



Vo 



tWo 



1 + eir-r-RJ/aR I _|_ ^{r-ri)/a, 



(66) 



with r/j = 2r^(16)3 and rj — 2rj(16)3. A uniform 
charge density with a radius r^ = 2. x 1.3x (16)3 was used. 
The fit parameters were Vq = 150.4 MeV, r% = 0.784 fm, 
an = 0.897 fm, Wq = 44.4 MeV, r? = 1.005 fm, and 
aj — 0.745 fm. By calculating with and without the 
Coulomb interaction, a Coulomb correction can be ob- 
tained by the method shown in Eq. [39l Alternatively, 
we can make the correction with the method Faldt and 



Pilkuhn 43[ by modifying the integral over the profile 
function. As seen in Fig. [5] the results are nearly identi- 
cal over most of the angular range. 
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FIG. 23: The ratio of the ''°Ca-'*°Ca cross section to the 
Rutherford cross section at 3.387 GeV using the appropriate 
parameters from Ref. [SO ] for this energy. 



The result of the ^^0-'*°Ca scattering is shown in Fig . 
[22]and compared with data of Roussel-Chomaz et al. [7l| . 
Again one sees the importance of the center-of-mass and 
Coulomb corrections. The agreement at forward angles 
is moderately good but is considerably worse at larger 
angles. For another treatment of this reaction at this 
energy see Ref. |lll| . 



J. ''°Ca-*"Ca scattering 

Figure [23] shows the cross section for '^^Ca-'^'^Ca scat- 
tering at a kinetic energy of 3.387 GeV. There is no data 
to compare with and if there were they would be dom- 
inated by pure Coulomb. In this case there are 1600 
factors in the product in Eq. [T] which, if expanded in 



40/ 



Ca- Ca scattering 
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multiple scattering as was done in the paper by Franco 
and Yin [6] would give 2^^°° « 4 x lO-^^^ terms. The 
calculation is shown to illustrate the fact that, even for 
these medium-heavy nuclei, the center-of-mass correction 
remains important. We see that the effect of the CMC is 
not very different from the carbon-carbon scattering and 
is as large as a factor of 5. Since the basic CMC effect 
is the order of 1/40 in this case it is important to under- 
stand how many times it enters in to the calculation. For 
b small the profile function is unity and only a relatively 
few factors are needed to achieve this result. 



K. 



^He-^^C scattering 
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"He+'^C 41.6 MeV/A data renormalized (0.9) 
Pudliner density (R =1.84 fm) (Exp 1.912 fm) 
Dash: With Coulomb, "no Shell Model 
Dotted: No Coulomb, no Shell Model 
Solid: With Coulomb and SM correlations 
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He + '^C 38.3 MeV/A data renormalized (.85) 
- Pudliner density (R =1.84 fm) (Exp 1.912 fm) 
Dash: With Coulomb, ""no Shell Model 
Dotted: No Coulomb, no Shell Model 
Solid: With Coulomb and SM correlations 




FIG. 25: ^Re-^'^C scatteringat 41.6 MeV/A using the ^He 
density from Pudliner et al. 0|. The data are from Al-Khalili 
et al. \M 



where 

e„ = 5.72753 fm^^ and /„ = 1.65566 fm"^ (70) 

A "shell model" correlation was (optionally) included by 
coupling two i = I neutrons to zero total angular mo- 
mentum which results in a factor in the density of 



2(' 



cos^ 9r) 



(71) 



FIG. 24: ^He-^2Q scatteringat 38.3 MeV/A using the "^He 
density from Pudliner et al. [J| . The data are from Lapoux et 
al. [6J| . The expe rimental determination of the proton radius 
is from Ref. [112( |. 

The Pudliner density for ^He Q was represented by two 
auxiliary functions (one for protons and one for neutrons 
following Eqs. [21] and [22] The proton auxiliary density 
was given by 



Vp(r) = 0.273776 



ri5&i6e- 



15! 



0.726224 



r°Cpe 



8! 



with 



(67) 



where Oj^ is the relative angle between the two valence 
neutrons. This factor was carried as a weight in the 
Monte Carlo calculation so the calculations with and 
without the factor can be done in the same run. Figures 
[Ml and [^ show the results for the calculation of ^He-^^C 
scattering. It is seen that the effect of the shell-model 
correlation is very small. The Coulomb correction plays 
a significant role even for these low-Z nuclei. Perhaps a 
part of the reason for this is that the correlation is only 
among two of the six particles. 



IV. DISCUSSION 



bp = 14.3466 fm"^ and Cp = 4.7125 fm"^ 
and the neutron auxiliary density by 



The Glauber approximation can fail for a number of 
reasons not necessarily associated with the basic assump- 
tions. These reasons include: 



'-pnir) = 0.208969 



r"e„ei5e- 



14! 



0.791031- 



'fne 



-Ur 



2! 



(69) 



• At low energies the Fermi motion may cause signif- 
icant corrections to the fixed-nucleon approxima- 
tion. 
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• The double spin flip in the nucleon-nucleon inter- 
action which has been neglected may require a sig- 
nificant correction if the isospin constraints are not 
sufficient to eliminate it. 

• The single spin flip (occurring an even number of 
times) may become important for sufficiently large 
angles. 

• At high energies the coherent production of mesons 
constitutes an additional inelastic channel which is 
beyond the Glauber approximation. 

• At small impact parameters there should be signifi- 
cant corrections to the nucleon-nucleon interactions 
because of the higher nuclear densities. 

• Correlations from the shell model are of the same 
range as those from the center of mass and may play 
a role as the nuclear penetration becomes greater. 
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FIG. 26: The factor of the kernel bJo{qb) which represents the 
weigliting of the nuclear profile function (top panel) and the 
nuclear profile function for a— Pb (lower panel). 

In looking over the comparisons with data presented 
in the paper it is seen that often there is a forward re- 
gion where the prediction is moderately good (sometimes 
quite good) followed by a more-or-less sudden transition 



to a poor agreement. Perhaps the best example of this 
is the scattering of a-particles from 208pj^ j-ggg pjg^ [T7| . 
The transition is seen to take place at 8° at 288 MeV, 
7° at 340 MeV, 5° at 480 MeV and 4.5° at 699 MeV. 
These angles and the corresponding momentum trans- 
fer and values of -t are shown in Table [TTl All of the 
transition points are seen to correspond to a momentum 
transfer of about 1 fm"! [-t=0.04 (GeV/c)^]. For the 
i6q_i2(^ scattering (see Fig. [20)1 the transition angle is 
about 4° which corresponds to a momentum transfer of 
1. 01 fm-i. For i^O-^^O scattering (Fig. ^ the transi- 
tion angle is around 5°, corresponding to a momentum 
transfer of 1.28 fm~^. For the ^^0-^"Ca scattering (Fig. 
[22)) the transition angle is about 2.5° which corresponds 
to a momentum transfer of 1.06 fm~^. The remainder of 
the cases studied are shown in Tabic |ll] as well and show 
a similar effect. 

Figure [26] shows that a.t qc — I fm^^ for the momen- 
tum transfer the interior of the nucleus is sampled. So 
qc — I fm~^ does not seem to correspond to the onset of 
the nuclear matter being probed. 
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Appendix A: Partial-wave projection 

We make the partial wave projection from the formula 

— - V(2£+ 1)(S'« - l)Pe(cos0) = ik / hdb,Uqh)G{h). 

(Al) 
where q^ = 2fc^(l — cos 9). In order to make this projec- 
tion it is useful to have the expansion 



Mqb) = Y,he{kb)Piicose). (A2) 

From Ref [ul 9.1.79 we have 

OO 

Jo {qb) = J^ {kh) + 2^Jl ikh) cos nO. (A3) 

n=l 

In order to get the desired expression we need the expan- 
sion 



cos n9 — y ^ an^iPijcosi 



(A4) 



e=o 
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Case 


TLab (GeV) [Ref 


] Oc{ deg) 


gc(fm-i) 


ql {GeV/cf yi- 


- rriT — mp(GcV) 


Figure 


a-2"«Pb 


0.288 [69] 


8 


1.04 


0.042 


0.28 


fTTl 


a-208pb 


0.340 [69] 


7 


0.99 


0.038 


0.33 


|1T| 


a-^o«Pb 


0.480 [69] 


5 


0.85 


0.028 


0.47 


m 


a-208pb 


0.699 [69] 


4.5 


0.93 


0.034 


0.65 


m 


16o_12c 


1.503 [71J 


4 


1.02 


0.041 


0.64 


|20| 


16q_16o 


1.120 [67, 68] 


5 


1.28 


0.064 


0.55 


|21| 


i^O-^Ca 


1.503 [71J 


2.5 


1.06 


0.044 


1.06 


m 


12(^_12c 


1.016 [70] 


4.5 


0.96 


0.036 


0.50 


|18| 


Q-1«0 


0.240 [36| 


7 


0.66 


0.017 


0.19 


m 


a-^«0 


0.400 [72] 


7 


0.89 


0.031 


0.32 


m 


12(^_12c 


1.449 [73] 


4 


1.01 


0.040 


0.71 


|19| 


12c_12(. 


1.620 [74| 


3 


0.80 


0.025 


0.80 


m 


12(^_12^ 


2.400 [73] 


2 


0.66 


0.017 


1.17 


[H 


a-«Ca 


1.37 [75] 


8 


2.18 


0.184 


1.23 


|16| 


a-i^C 


4.20 [79, 80] 








2.87 


m 


a- C 


1.37 [76] 








0.99 


[13] 


a- a 


2.55 [77] 








1.18 


El 


a- a 


4.20 [78] 








1.87 


El 



TABLE II: Critical values of the angles and their equivalents in momentum transfer. Also given is the energy in the center of 
mass available for coherent meson production. The last four cases have no critical values listed because all data points available 
have a momentum transfer greater > 1 fm~^. 

The ^He-^'^C cases have not been included since it is not possible to claim agreement in the forward direction without absolute 

data. 



which actually cuts off at £ = n as we shall see. With the 
use of DeMoivre's theorem on can see 



cos 



m=0 \ / 



E 



m^even 



;i™ + (-i)™]sin™6'cos"-™6' 

(A5) 

(-1)^ sin™ 61 cos"-™ 6* 



where is the binomial coefficient. With k = to/2 

\ n 



2k 



fc=0 V " 



2fc 



2 a\k „„„n-2k , 



fc=0 



(-l)'^(l-cos^6')'^'cos 



^^ / 2fc \ (_i)fe(_i), ( J ] cos"-2'=+2j" 



k 



With 771 = n — 2k + 2j so that ni = n, n — 2, n — 4, . 
71 — m — 2j 



Er^"' ' (-1)^ 



n — m 



J 



= Y^ Cn,m COS™ 0. 

Thus we see that cos 7i9 can be expanded in powers of 
cos 9 with powers of the same parity as n with maximum 
power n. We can now use the int egra l of a Legendre 
polynomial over a power of x (Ref. |ll3J | 8.14.15) 

/ N m , _ 7r3 2-"T(l+TO) 

^^Pe{x)x dx ^^^ ^^^_ ^^^^^^^ +lm+l)- '^™' 

(A6) 
Note that this integral is zero if £ and to, are not of the 
same parity or if £ > to,. 

We can write the expression for the a„_^ as 



1, 



£<m<7i,m~7i even 



E 



(A7) 



This expression was coded and the coefficients calcu- 
lated in this way; it works but the terms for different 
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powers of cos 9 get very large and cancel in the sum so, 
even with double precision, one is limited to n of the or- 
der of 35 for small values of £ which is not enough for our 
problem. The values ioT n = £ and other values of £ close 
to n are always calculated correctly. A useful check is 



y^ an,i = 1 



i=0 



(A8) 



which follows from Eq. IA4I with 6 = 0. 

However, there is a much easier way to get the coeffi- 
cients needed. Using the trigonometric identity 

cosn6' = 2cos(n - 1)6* cos 6* - cos(n - 2)0 (A9) 

and the recursion relation for the Legendre polynomials 



2^ -Hi 



we find 



r 

/ cos nOPi (cos 9) sin 9d9 = bn,£ 
Jo 

= ^7^ [^&n-l,£-l + (^ + l)bn-l,i+l] - bn-2,t- (All) 

Since a„^£ = ^^Y^bn^e we have 

2£an-i,i-i , 2{£ + l)an-i,e+i /Air,\ 

««,£ = — —n — ; 1 ^„ , „ an-2,e (A12j 



2£ - 1 2£ -t- 3 

with the conditions 



and 



o-ofi — 1; ai,i — 1; 



a-n.e = 



(A13) 



(A14) 



if n and £ are not of the same parity or if ^ > n. All values 
of these coefficients can be quickly calculated using the 
recursion relation IA12I without any numerical problem. 
Using these values then 



ho{kb) = Jl{kh) + 2 ^ JUkb)a„ 



hi{kb) 



2Y, Jn{kb)anj; £>0 



(A15) 



(A16) 



Jo{x) and Ji(x) are calculated to one part in 10^ and the 
recursion relation for Bessel functions is used for fixed b 
to obtain the higher values of n. Just beyond n = bk, the 



magnitude of Jn{x) decreases rapidly and the recursion 
relation fails shortly thereafter. The recursion is cut off 
when Jn{x) is less than 10~^. 

The functions he{kb) have some interesting and poten- 
tially useful properties. For example starting from Eq. 

2] with 



Mqb)Mq'b)bdb = d<^^Hq-q') 



t^ = 2S{q^ q'^) = Isix - x') (A17) 



one can show that 



2£-|- 1 
hg{v)ht{y)ydy= — - — Siji (A18) 



Once the functions hi(bk) have been found, the pro- 
cess can be reversed, i.e. we can find the nuclear profile 
function for any given set of matrix elements. First from 
the inverse Hankel transform we have: 

G{b) = ^ r F(q)Uqb)qdq (Al9) 

*^' Jo 

Since the amplitude is assumed to fall off rapidly we 
can replace the upper limit with a finite, but large, one, 
namely 2k so that 

G{b) « 1 / Fiq)Mqb)qdq (A20) 

Changing the integration variable to u = q^ we have 

1 f^''' 
G(5)«— y FiV^)Joiqb)du (A21) 

Now changing the integration variable to a; = cos^ we 
have 

G(6)«— y Aix)Mqix)b)dx (A22) 

where A{x) — F[q{x)] and 

A{x) = —J2(2£+l){Se-l)Pe{x) (A23) 

e=o 

Using Eq. IA2l we have 



G{b)^-J2(Se-l)he{kb) 



(A24) 



e=o 
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Appendix B: Variational density for *He 

One can solve for the ^He wave function by Monte 
Carlo Green's Function methods @) Q • ^^ this case the 
walkers either represent a probability proportional to the 
wave function itself (hence the walkers do not give the 
density needed for our scattering problem) or, if one uses 
importance sampling then the walkers represent the prod- 
uct of the trial wave function and the true wave function 
and thus is more suitable for representing the density. In 
the work here we will use a variational trial wave func- 
tion where the walkers represent directly the density to 
calculate a — a scattering. 

Varga et al. [3J| performed Monte Carlo calculations 
using a Metropolis sampling of a ''He wave function ob- 
tained with Green's Function Monte Carlo Methods. We 
will make a comparison between the method given previ- 
ously and a simplified version of this type of calculation 
based on the variational algorithm. 

The variational method operates with the estimator of 
the energy given by 



Et — 



/dRV;|,(R)(T-HF)V>T(R) 
/dRV'?.(R)V'T(R) 

JdRpiK) 



(Bl) 



where R represents the set of coordinates which describe 
the system. Using the Metropolis algorithm to repre- 
sent the density, it is automatically normalized. The trial 
wave function is assumed to depend on some number of 
parameters and, upon varying these parameters the min- 
imum energy achieved is guaranteed to be greater than 
or equal to the true ground state energy of the system. 
The trial wave function normally is expected to give a 
good representation of the true wave function but there 
is no guarantee of that. 

In order to calculate with a realistic wave function we 
use a trial wave function which results from the follow- 
ing variational calculation for '^He. Since the true wave 
function must be translationally invariant it can depend 
only on the six relative coordinates, r^ = r^ — r^. 

We choose the form 



^(ri,r2,r3,r4) = J|/(?" 



(B2) 



j>i 



where 



\rij\ and the function, /(r) is arbitrary at 



this point. It will be chosen with a number of parameters 
to be selected to minimize the energy. We can evaluate 
the kinetic energy needed as 



i=l 



Vf-0(ri,r2,r3,r4) 

V'(ri,r2,r3,r4) 



-2E 



j>i 



fin,) , , fin,) 



finj) r^jfinj) 



E' 



Vtk 



f{n,)f'{nk) 

ri]rikf{rij)f{rik) 



(B3) 



The values of the three indices in the last sum are: [123], 
[124], [134], [213], [214], [234], [312], [314], [324], [412], 
[413] and [423]. They can be obtained with nested loops 

{i = 1,4[.? = 1,3; J ^ i{k = J + 1,4; fc ^ t)]} (B4) 

With such a form it is relatively easy to calculate the 
kinetic energy since one only needs to be able to calculate 
the function, /, and its first and second derivatives. The 
derivatives can be calculated numerically if need be. 
We took for the form of / 



/(r) = (l-e— )- 



(B5) 



The first factor vanishes at the origin and provides the 
effect of a (very mild) repulsive correlation. The rest of 
the function has the proper asymptotic form with the 
constant b avoiding the singularity at the origin. 



We used a modified version ISTJ the Malfliet-Tjon |ll4| 
potentials since we are neglecting spin in this calculation. 
The variational calculation over-binds *He by about 7 
MeV (35 MeV instead of 28 MeV) so we increased the 
strength of the repulsive part of the modified Malfliet- 
Tjon potentials (by 6%) in order to have a more realistic 
binding. 

The variational calculation for the energy was carried 
out using 100,000 walkers. It was found that the energy 
was not very sensitive to the parameter b so it was fixed 
at 1.0 fm. Then, for various fixed values of c, the energy 
was calculated as a function of a. For the fixed values of 
b and c the rms radius is a unique function of a. We used 
a = 0.203 fm~^ and c — 1.0 fm~^. The rms radius is 
a crucial parameter in the scattering and we obtained a 
value of 1.44 fm. The calculation was repeated 10 times 
to give a collection of one million walkers to use in the 
scattering calculations involving the a particle. 
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The coordinates needed for the scattering calculation 
were taken from the Metropolis algorithm every 2000 
steps to be certain that each realization was independent. 
One million realizations of the helium density were calcu- 



lated and written out to files. When the scattering was 
calculated, the coordinates of each nucleus were drawn 
from a nuclear representation chosen randomly from the 
pool of one million possible nuclear realizations. 
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